Bath Assisted Cooling of Spins 
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A suitable sequence of sharp pulses applied to a spin coupled to a bosonic bath can cool its state, 
i.e., increase its polarization or ground state occupation probability. Starting from an unpolarized 
state of the spin in equilibrium with the bath, one can reach very low temperatures or sizeable 
polarizations within a time shorter than the decoherence time. Both the bath and external fields 
are necessary for the effect which comes from the backreaction of the spin on the bath. This method 
can be applied to cool at once a disordered ensemble of spins. Since the bath is crucial for this 
mechanism, the cooling limits are set by the strength of its interaction with the spin(s). 

PACS numbers: 05.30.-d,76.20.+q 



Cooling, i.e. obtaining relatively pure states from 
mixed ones, is of central importance in fields dealing with 
quantum features of matter. Laser cooling of motional 
states of atoms is nowadays a known achievement l^. 
The related problem of cooling spins is equally known: 
it originated as an attempt to improve the sensitivity 
of NMR/ESR spectroscopy 0, H H S H 0, since in 
experiments the signal strength is proportional to the 
polarization. Recently it got renewed attention due to 
realizations of setups for quantum computers The 
very problem arises since the most direct methods of 
cooling spins, such as lowering the temperature of the 
whole sample or applying strong dc fields, are not fea- 
sible or not desirable, e.g. in biological applications of 
NMR. Indeed, at temperature T — IK and magnetic field 
B = IT the equilibrium polarization of a proton is only 



tanh 



hfj.B 

2fcBT 



10 since the ratio /i 



frequency 
field 



is equal 



to 42 MHz/T. For an electron /i is 10^ times larger and 
for ^^N it is 10 times smaller. The weak polarization can 
be often compensated by a large number of spins, but 
for some NMR-isotopes the natural abundance is too low 
(0.36% for i^N). 

Over the years, several methods were proposed to at- 
tack the problem of small polarizations. The polarization 
is generally increased via a dynamical process and it is 
used before relaxing back to equilibrium H, H, IE IM ^ Q - 
Specially known are methods where a relatively high po- 
larization is transferred from one place to another, e.g. 
from electronic to nuclear spins |M Q IS - In this 
respect electronic spins play the same role as the zero- 
temperature bath of vacuum modes employed for laser 
cooling of atoms [iJl (this bath is typically inadequate 
for cooling nuclear spins, but can be employed to study 
cooling of atomic few-level systems in the context of opti- 
mal control theory |^). Polarization transfer was studied 
in various settings both theoretically and experimentally 
1^ IE ^ IE IS' However, this scheme is limited — besides 
requiring an already existing high polarization — by the 
availability and efficiency of the transfer interaction. A 



related method, polarization compression, consists in ma- 
nipulating a set of n spins in such a way that the polar- 
ization of one spin is increased at the expense of decreas- 
ing the polarization of the remaining n — 1 spins. These 
spoiled ones can be recycled and used again 0. Since 
spins are cooled one by one, a long time and carefully 
designed inter-spin interactions are needed for cooling a 
large ensemble. 

Here we propose a mechanism of cooling which only 
uses the most standard setting of NMR or ESR physics 
[EIE0|: spins- ^ under the action of external field pulses 
coupled to a thermal bath at the same temperature. 
The bath is needed because external fields alone cannot 
achieve cooling ^El- However, we assume neither that 
the bath is under any direct control, nor special con- 
straints on the bath-spin interaction: it is the standard 
one, widely studied in the context of decoherence. We 
show that, rather than being a hindrance in quantum 
system manipulations, the bath is capable of producing 
ordered effects on the spin, which can cool it down to 
very low temperatures (~ 1/xK for a proton) in a finite 
time. Two factors are crucial: the backreaction of the 
spin on the bath and the generation of transversal com- 
ponents (coherences) during the cooling process. Since 
the effect is generated via the bath, one can cool at once 
a completely disordered ensemble of spins. 

The model we study is well known 

dmilEIIElIi: a 

spin-i with energy levels ±^Hil couples to a bath, mod- 
eled by a set of harmonic oscillators with creation and 
annihilation operators and dk- The total Hamiltonian 
reads 
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[ai,al] 



6ki. (1) 



Here ujk are the bath frequencies, a^^y^z the Pauli opera- 
tors, and X — J2k9k{^k + the collective coordinate 
of the bath. The interaction is chosen assuming that 
the Ti-time, connected to relaxation of the average (az) 
is very large (infinite) 0, IE 0- 

The gk are couplings 
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parametrized via the spectral density function J(lo): 

J{^)^Y./k5{^-^k). (2) 

In the thermodynamic hmit the bath modes are dense 
and J{uj) becomes a smooth function determined by the 
physics of the system-bath interaction "TTl. The oscilla- 
tors can represent real phonons or stand for an effective 
description of a rather general class of thermal baths ■ 
Let us recall how the model is solved |l2(: az is 
conserved, while ak{t) = e-*"'=*afc(0) -I- (5"*'^'=* - l). 
This leads along with Eqs. (PEJ to 



X{t)^f^{t)-dzF{t), (3) 
m EE V afc[4(0)e-'^-* -|-afc(0)e— (4) 



J[uj){t- 



sin Lot , 



(5) 



where f}{t) is the quantum noise operator, and where 
Fit) = ^F{t) quantifies the backreaction of the spin on 
the collective operator of the bath. This effect, not rele- 
vant for decoherence as such, is crucial for our purposes. 

We assume that at the initial time t = the common 
density matrix of the bath and the spin is factorized: 



PiO) 



-l3Ho 



tr e-/5^o 



Ho 



^^JujJka\ak, (6) 



where T = is the common temperature {k-Q = 1). 

p(0) describes the spin prepared independently from the 

bath and then brought in contact with it at i = 0, e.g., 

by injection of the spin into a quantum dot or by creation 

of an exciton through external radiation. 

As follows from Eq. (O, fi{t) is a Gaussian operator 

< 

with {fi{t)) = and time-ordered correlator {fj{t)f]{0)) = 
^{t) — iF{t), t > 0, where (...) is taken over the initial 
state (|nj, and where 



^ , , 1 — COSLUt , hlO , , 

dujJ{uj) 2 coth— . (7) 



The Heisenberg equation of the spin, ha± — i[H,a±\ 
with a± — cTx^i(Ty, aza± = ±(T±, is solved as 



<7±(<) 



<7±(io), (8) 



where V is the time-ordered exponent. Defining £tA — 
^itH/h^^~itH/h Qj^g derives 

xitiM.t) = F{t2) - F{ti) + F{ti +t)- F{t2 + t), (10) 

Eq. is the standard formula for the average of 

a Gaussian operator. The factor e^'''^*' leads to de- 
coherence 12], since due to Eqs. (|H1 EJ, {cr±{t)) — 



e-5(*)±*f^*(i7_j_(0)) for a general factorized initial state. In 
this simplest situation the backreaction factor F, prop- 
erly obtained already in ^3>Q|> cancels out. In general, 
F can shift the spin's frequency O as seen below. 

The action of external fields on the spin amounts to 
a time-dependent Hamiltonian H{t) = H + h{t)(j. In 
the pulsed regime 0, S S E H 13 M<) difi'ers from 
zero only for very short intervals of time 6 being there 
very large, h{t)S ~ 1, to achieve a finite effect. As a 
consequence, terms oc CTz in can be neglected during 
the time-interval 6. A single pulse can perform an ar- 
bitrary unitary transformation in the space of the spin 
(rotation of the Bloch vector (it)). We parametrize it as 
U = e'«'^(*)^A, (0 < (/), V' < 27r, < 61 < 7r/2) 



U 



-i4> 



id -e-*'^sini9 



e"^ cos 



, VA=UAU\ (12) 



As a first example, we take ohmic interaction l.llj: 

J(tj) =7cje-'^/^. (13) 

where 7 is a dimensionless coupling constant, and where 
r (usually ^ f2) is the bath's response frequency. Eqs. ((SJ 



Elinj imply ^{t) =7 In 



r^(i+e) vi+r-'t'^ 



r(i+e-ie r t)r(i+e+ie r t) 



F{t) = 7 [n - arctan(ri)] , 9 = 



kBT 



(14) 



where T is Euler's function and Q is the dimension- 
less temperature. The time-scale l/7r of backreaction 
F{t) is temperature-independent as opposed to the deco- 
herence time. For low temperatures 9 <C 1: e^'»^*^ — 
(1 -I- t2r2)-7/2^ while for 9 > 1, e'^^*) starts as agaus- 
sian, but continues as e"*/^^ with 7^ = h/{2^T) [l2| . 

Cooling amounts to make the final polarization (a'z)f 
more negative than the initial one, (ffz)! = — tanh ^(3h^l 
(Eq. 10) implies {ax)i = {^y)i = 0). A single pulse cannot 
achieve cooling since it sees the initial local equilibrium 
state of the spin, and then according to the no-cooling 
principle \u\ it can only heat the spin's state up: for an 
arbitrary pulse Vi applied at time i, {crz)i = {^t Vi Oz) — 
((72)i cos2?9i > {(Jz)i (recall {az)i < 0). Thus we have 
to employ at least two pulses. The final polarization 
after one pulse at t and one at i -I- t, P = |((5'z)f| — 
\{£tVi£rV2^z)\, reads from Eqs. 

(o'z)t = ((Tz)! cos2i3i cos2??2 S2 sin2i3i sin2z?2, (15) 
S2 = -e-«M3?{e'""+^"^ ((a,)iCosx + isinx)}, (16) 

where x = xlOj^i''') defined in Eq. I|l()|l . and 

ct2 — V'l — '^2 — 4>i — <t>2 arises from Eq. \i'2.l . There are 
now two factors that come from the bath: e^^^"^-* in S2 
accounts for the decoherence in the time-interval (t, t-\-T) 
of transversal terms generated by the first pulse, while x 
is the backreaction factor from Eqs. IpH 110(1 . 

Though the finite-i situation can be of its own interest, 
for all results below we set ^ 1 (a mild condition, since 
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TABLE L Maximum value of | {(Jz) \ for 2, 3 and 50 x 2 pulses 
and various values of 7 and Q = k-BT/ifiV) (for the ohmic 
case) . The lowest two lines represent the final collective polar- 
ization \m{ \ of a strongly disordered ensemble with fio = and 
arbitrary large disorder dispersion d under spin-echo setup of 
2 pulses plus a tt pulse and a sequence of 50 such sets of pulses. 



FIG. 1: Final polarization P = l(o"z)f | given by Eqs. II15II16II 
versus dimensionless time between two optimal pulses. Left: 
ohmic interaction. Bold curve: 7 = 5, O = 0.2. Normal: 7 = 
2, O = 0.5. Dashed: 7 = 1, O = 1. Right: 1// interaction, 
"if > 10^ and Of/^f = 0.01, 0.1, 1 (bold, normal, dashed). 
In both figures ((Tz)i and are neglegible. 
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1/r is typically the shortest time-scale), since this makes 
the outcome independent on the details of the initial state 
preparation. In this ergodic limit, the initial condition 
p(0) oc e~^^^ defined by Eq. © is equivalent to the 
overall equilibrium preparation /5cq(0) cx e^^^ |l9l |. 

In Eq. (|16|l . S2 can always be made negative by tuning 
a2. Minimizing ((T2)f over '(?2 produces min [((T^)!, S2]. 
If the initial polarization is already high, |((T2)i| > |s2|; no 
pulses should be applied, since they only heat the spin 
up. However, in the relevant situation {cTz)i — 0, the 
minimum {az)i = S2 is reached for — §2 ~ ^ ~ ■j- 
Altogether, using Eq. ((Tni) and ^ F, yields 

{az)f = S2 = -e"^(^^ sin [7 arctan(Tr)] . (17) 

The choice of optimal pulses, which have to be coherence 
generating, can be a |-pulse along the x-axis followed by 
a — ^-pulse along the y-axis: ViUz^x = '^y,x, T^i^y = ~^z, 
and V20'z,y = ^x,y, 'P2<^x = —<^z- Fig- n E^nd Table I 
show that max,-|((Tz)f I can approach its maximal value 1. 
Transversal components generated by the two pulses will 
decay after a time T2, and the spin will be described by 
a Gibbsian at a temperature lower than the initial T . 

The origin of this cooling effect is in shifting the spin's 
frequency by the factors F ~ 7r and Xi recall Eqs. 
I10|l . which arise from the (via the pulses) enhanced back- 
reaction of the spin on the collective coordinate X of the 
bath. The generation of coherences by the first ^-pulse 
is necessary to couple ct^ to the bath, which so to say 
"thermalizes" {a^) under the shifted frequency. Cooling 
is achieved via the proper pulses, still it decreases with 7 
(weaker backreaction) and with (3 (larger decoherence). 
The time r between two pulses can be neither too short 
(x is visible on the time-scale l/(7r)), nor too long, since 
otherwise decoherence will diminish the influence of the 
first pulse. As shown in Table I, the cooling improves 
by i) applying three successive pulses and optimizing over 
their parameters; ii) applying two pulses, waiting for a 
time > 72, so that the transversal components decay, 
((7±) 0, applying another two pulses, and so on n 
times. The final K^z)! is maximized over all free param- 
eters. It appears that this numerical maximization can 



be done locally, i.e., by maximizing the output \{crz)\ af- 
ter each pair of pulses. This "greedy" optimization shows 
up also in cooling via a zero-temperature bath 0. 

Inhomogeneous broadening. Many experiments in 
NMR/ESR are not done with a single spin, but with 
an ensemble of non-interacting spins having random fre- 
quencies f2 due to action of their environment or due 
to inhomogeneous external field 0, 0- The collective 
variables are obtained by averaging the corresponding 
expressions for a single spin: rrif = ^ dQ.P{Vl){az)i. As- 
sume that the distribution of U, is gaussian with average 
r^o and dispersion d: P(il) cx e^^^"^") Averaging 
over P(ri) the term e*^'^ in Eq. (|16|l produces a factor 
^ ^~dT'/2^ strong decay on times 7^* oc \/y/d. After 
this decay, S2 ^ and any two pulses will only heat the 
ensemble up as seen from Eq. (|15|) . 

It is however possible to employ the spin-echo phe- 
nomenon and cool, i.e., increase the collective final po- 
larization |mf I as compared to the initial even for a 
completely disordered ensemble with Tj* being very short: 
Apply precisely in the middle of the two pulses an addi- 
tional TT-pulse in x-direction: V-nC^y^z — ~^y.z, V-nCfx — 
dx, and work out TOf = ^ dD. P{Vl){Et'PiEr'PTrEr'P2(yz)'- 

TOf = — TOi COS 2-di cos 2^2 + S3 sin 2-di sin 2i?2, (18) 
S3 = -e-''«(^)+«(2^)5R{e^"^ (micosxa - ^sinxs)} , (19) 
X3 =x(0,T,t)-x(T,2T,t), (20) 

where 03 = — ipi — (f>2 — tp2- As compared to (|lt)|l . 
both the decoherence e~45('^)+?(2-r) ^^^^ ^j^g backreaction 
X3 term are different. In the gaussian regime ^ oc 
decoherence is absent, while the exponential regime ^ oc t 
is unchanged, e-4«(^)+«(2r) ^ g-^^^r) Due to 

the Vtt pulse, the 7^*-decay has been eliminated, no term 
like e~*^* in Eq. (|16ll appears here. Now flo and d enter 
only via rui. The structure of Eqs. H18I20|I is close to 
the one of Eqs. (|15I16|I . and the optimization over 'di,'&2 
goes in the same way. To facilitate comparison, we take 
r^o = 0, thus TOi = 0, and disorder strength d arbitrary 
large, recalling in addition Eqs. (|10II13|) and the ergodic 
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condition tT ^ 1, we get 

\mf \ = e"^«(^)+«(2^) sin{7 [2arctan(Tr) - arctan(2Tr)] }, 

where we already inserted the optimal values iSi — "§2 ^ 
— ^ = j- The choice of optimal pulses can be the same 
as for the two-pulse scenario. The maximal |77if| can ex- 
ceed 0.4 for sufhciently strong coupling and/or low tem- 
peratures. The results improve by applying a sequence 
of three (spin-echo) pulses separated from each other by 
a time much larger than 72, see Table I . 

A 1// spectrum is another relevant situation of the 
bath-spin interaction recently observed in a two-level sys- 
tem (spin) of charge states in Josephson-junction cir- 
cuit (Cooper-pair box) The spin's interaction with 
the bath of background charges is modelled via 0, 0| 
Jf{ijj) — e^^l^ 9{uj — A), where &/ is the coupling con- 
stant, 9 is the step function, and where T and A are, re- 
spectively, the largest and smallest frequencies of bath's 
response. The upper frequency F is not relevant: the in- 
tegrals in Eqs. |7|) converge for F ^ oo |0. Within 
the realization of J^J: F ^ cx) (» 10"Hz), A = 314Hz, 
Ml = 34/ieV, while the relevant dimensionless coupling 
constant is very large 7/ = hf/h? = 2.3 x 10"'^^. Thus 
the cooling effect will be especially visible here. The 
spin-boson model does not describe all the aspects of 
the decoherence in Cooper-pair box; see [T^ for more 
elaborated approaches. Still some important features 
are reproduced adequately 0, 0|, and this motivates 
us to work out Eqs. [T)) and apply them to the 
two-pulse cooling situation described by Eqs. IjlSI 
under the following conditions, i) high temperatures 
0/ — k^T/lhh) ^ 1 (pessimistic case!); ii) tA <C 1 
and 7/ 3> 1 (experimentally relevant regimes [13); ill) 
{az)i — due to large temperatures, iv) ftr ^ 7/A due 
to large 7/. In analogy to Eq. (fTTjl . the result is {a2 = f ): 

{az)f — e~^fy'^^"'f sin — |^ j , where y = 7/At is a di- 
mensionless time, and where we omitted terms 0[y^ /j'j]. 
As seen in Fig. ^ the polarization can increase from its 
original value zero to 0.997 for S//7/ — 0.01, which in 
numbers of Ref. corresponds to T = 15K. 

In conclusion, we described a new method for cool- 
ing spins due to common action of external fields and a 
bosonic bath, starting from the overall equilibrium state. 
The fields alone cannot cool ^3 j while the bath alone can 
generate only the standard decoherence ll^. As com- 
pared to existing methods [IjlESIEIEHUjthe present 
one assumes neither already existing high polarization 

IE Q nor controlled spin-spin or bath-spin in- 

teractions U, nor a low-temperature bath It works 
even for very weak dc fields and applies to an ensemble 
of spins having completely random frequencies (strong 
inhomogeneous broadening). We are not aware of other 
methods achieving such a goal. The spins are cooled 
at once (not one by one) and the cooling process takes 
a time shorter than 72. Together with the overall effi- 



ciency of the method, see the figures and the Table, these 
features are encouraging for applications, e.g., in NMR 
spectroscopy. Our basic assumptions are a decoherence 
time T2 much smaller than the energy relaxation time 
7i , and the availability of sharp and strong pulses acting 
on the spin. A long 7i time characterizes other methods 
0, 1^ , while strong and short pulses were used for a clean 
demonstration of the effect, which probably survives for 
other types of pulses. 

The origin of the present mechanism lies in shifting the 
spin's frequency due to backreaction of the spin on the 
bath. This dynamical effect requires a non-perturbative 
treatment of the bath-spin interaction and is usually 
missed by markovian approaches |18|. It operates on a 
specific time-scale and allows to cool the spin provided 
the proper, coherence generating, sequence of external 
pulses is chosen. The cooling is efhcient already for small- 
to-moderate bath-spin couplings, and is especially visible 
for situations where a strong bath-spin coupling is inher- 
ent (l//-noise). In this experimentally realized situation 
[lij the cooling mechanism is expected to be feasible. 

We thank P.A. Bushcv for interesting discussions on 
cooling. This work was supported by FOM/NWO. 
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